Abstract. The isolated neutron star candidate RBS 1223 
was observed with the Advanced CCD Imaging Spectrom- 
eter aboard the Chandra X-ray observatory on 2000 June 
24. A timing analysis of the data yielded a periodic mod- 
ulation with a period P = 5.1571696+^* ^ s. Using 
ROSAT HRI archived observations we detected a period 
P = 5.1561274 ± 4.4 x 10~ 4 s and determined period 
derivative P = (0.7-2.0) x lO" 11 s s" 1 . The detection of 
this period and P indicates that RBS 1223 has a "char- 
acteristic" age of 6000 — 12000 years and huge magnetic 
field at the surface (B dlpole w (1.7-3.2) x 10 +14 G) typical 
[ for anomalous X-ray pulsars (AXPs). 
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1. Introduction 

Neutron stars (NS) supposedly represent a non— negligible 
fraction of ~ 1% of all stars in the Galaxy, with a total 
numb er as high as 10 s — 10 9 (e.g. Narayan & Ostriker 
199C). Isolated NS represent the bulk of the population 



and soon after their birth may appear as radio pulsars. 
This pulsar phase could last some 10 7 yr and these NS 
are observable, if the beaming conditions are favorable. 
At formation NS are very hot (T ~ 10 11 K), and they may 
give off thermal radiation while cooling. 

Depending on some model assumptions, a few hundred 
to several thousand accreting isolated NSs were expected 
to be detec table in the ROS AT A ll-Sky Survey (Treves & 
Colpi |1991|, Blaes & Madau p9§ Neuhauser & Triimper 



ct 
to 



199£). However, to date only 6-7 good candidates (Treves 
et al. |2000| , Motch |2000| ) for isolated NSs were discov 



ered in ROSAT data (RX J18 56.6-3 754, Walter |1996 
RX J0720.4-3 125, Haberl et al ~ 
Haberl et al. f998l 1RXS J130848.6 



1997| , RXJ0806.4-4123, 
212708= RBS1223, 



alogue, Voges et al. 1999) and studied as a candidate 
isolated neutron star (Schwope et al. 1999[ ) during the 
program of optical identifications of bright(> 0.2 PSPC 
cts s _1 ), high-galactic-latitude (\b\ > 30°) X-r ay so urces 
(ROSAT Bright Survey - RBS, Schwope et al. EOOOP . 
A pointed ROSAT HRI observation (Schwope 
al. 1999[ ) determined the position of the source 
be a(J2000) = 13 ft 08 m 48.17 s , <5(J2000) 
21°27'P7.5"(±1.6"). The corresponding Galactic coordi- 
nates are I = 339°, b = +83°. Deep optical observations 
with Keck-II failed to detect an optical counterpart down 
to a limiting magnitude of tub ~ 26, thereby placing a 
lower limit of ~ 10000 on the X-ray to optical flux ratio 
(Schwope et al. [l999[ ). 

In this paper we report the detection of pulsations from 
the isolated neutron star RBS 1223 based on recent obser- 
vations with the Chandra X-ray observatory and ROSAT 
archival observations. 



Schw ope et al. |l999j , RXJ1605. 3+3249, Motch & Haberl 2. Observations and Data reduction 



1999 and RXJ0420. 0-5022, Haberl et al. 1999). They can 



be exclusively described by their similar X-ray properties, 
i.e. they have soft X-ray spectra which are well repre- 
sented by pure blackbody emission with temperatures kT 
between 50 and 120 eV, and extremely high X-ray to opti- 
cal fl ux rat ios, log ( f x / f o v t) of 4— 5.5 (see tables in Treves 
et al. [2000j Motch |2000| ). 

Most probably, these best candidates constitute the 
bright end of the log N— log S distribution of X-ray de- 
tected isolated NSs (Neuhauser & Triimper 1999), sug- 



gesting that the smaller number of observed accreting old 
neutron stars than those theoretically expected is mainly 
caused by their larger space velocities compared to the 
previously assumed velocity distribution. This conclusion 
is supported by the recent d iscove ry of a proper motion 
for RXJ1856.6-3754 (Walter |200l| ). Two of these objects 
(RXJ0720.4-3125, RXJ0420. 0-5022) exhibit pulsations 
with remarkably long periods of 8.39 s and 22.7 s. 

The point source RBS 1223 (1RXS J130848.6+212708) 
was detected in the ROSAT All-Sky Survey (1RXS cat- 
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The point source RBS 1223 was observed with the Ad- 
vanced CCD Imaging Spectrometer (back illuminated S3 
chip) on board the Chandra X-ray observatory on 2000 
June 24, for a total of ~ 10 ksec. The source was placed at 
the farthest corner of the back-illuminated CCD S array, 
close to the read out of node 3 in order to mitigate pile- 
up and charge transfer inefficiency effects. Only a small 
part (1/8) of the CCD was readout to additionally avoid 
pile-up and increase the time resolution to ~ 0.5 s. 

We have used level 2 ACIS event data, retaining only 
'standard' grades 0,2,3,4 and 6 thus discarding events with 
very large pulse heights as due to cosmic rays. 

An inspection of an X-ray image in chip coordinates 
(rather than sky coordinates in which the dithering of 
space-craft has been accounted for) showed no 'bad' pixels 
or columns in the source data. Examination of the time 
dependence of the background did not reveal periods of 
high background during the observation. 

To identify point sources in the field (Fig. [I]), we used 
the routine CELLDETECT provided by the Chandra In- 
teractive Analysis of Observations (CIAO 2.0) software 
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Fig. 1. Chandra ACIS-S image of RBS 1223. Inner el- 
lipse corresponds to the 5 a radius of the source (available 
from source detection procedure, CIAO routine CELLDE- 
TECT) and elliptical annuli to the background covering 
an equal area. 



package. We determined the position of RBS 1223 to be 
a(J2000) = 13 /l 08 m 48.26 s ,(5(J2000) = 21°27'06.75", com- 
patible with the position found in previous ROSAT HRI 
observations ( see Sect. |]). 



2.1. Timing analysis 

2.1.1. Chandra X-ray observatory ACIS— S 

First of all, we performed a timing analysis of Chan- 
dra. X-ray observatory data in order to search for 
periodic modulations. We extracted ~ 7840 ACIS-S 
source photon events (Fig. Q), which we corrected for 
barycenter using the "axBary" programme, available at 



ftp : II asc.harvard.edu 



An applicati on of various methods (e.g. Scargle 1982 , 
Buccheri et al. 1983 ) for periodicity search indicates a 
highly significant detection of a periodic signal at 5.157 s. 

To evaluate the pulsation frequency more precisely, 
and find its uncertainty, we employed the Gregory-Loredo 
(GL) Bayesian method (Gregory & Loredo |l992j , |1993 , 
1996j ). This method is designed for the detection of a pe- 
riodic signal of unknown shape and period, assuming a 
Poisson sampling distribution. To compare different mod- 
els (how the data favor a periodic model of a given fre- 
quency / with m phase bins over the constant rate model) 
it calculates a frequency-dependent Bayesia n odd s ratios 
(O m (/)) (for details, see Gregory & Loredo 1992| ). 



We implemented the GL method, taking the maximum 
number of phase bins equal to 15 and surveyed the fre- 
quency range 0.01-0.25 Hz. The odds ratio as a function 
of frequency is shown in Fig. |2| (left panel). As a most 
probable value for the period we took the mode of the 
posterior probability distribution function (PDF). 

The maximum value, 0™ ax = 6.16 x 10 10 , is at / = 
0.19390681 Hz. We computed the uncertainty of / at 
68% confidence level ("posterior bubble") using the full 
Bayesian posterior PDF. It is equal to 5f =^'1^ /iHz. 
Thus, we finally derive the period of the detected pulsa- 
tions and the uncertainties to be 



P< 



5.1571696 



+ 1.57xl(T 



Chandra — o.±oi iuju_| 36xl0~ 4 

at the epoch of 51719.949676196 MJD. 

The light curve extracted at / = /o (Fig. || left panel) 
reveals one broad pulse per period. The semi-amplitude 
modulation, 

. C R max C Rrain 

^Chandra — 7775 T r< p ' 

where CR is the count rate, is equal to 23% ± 5%. As- 
suming that the detected signal is sinusoidal we estimate 
that the pulse d fraction intrinsic to the source is 20% ±2% 
(Brazier |1994[). 



2.1.2. ROSAT HRI 

The field of RBS 1223 has also been observed with the 
ROSAT HRI on June 13, 1997 (ROR 703348) and on Jan- 
uary 10, 1998 (ROR 704082) for a net exposure time of 
2218 s and 4938 s, correspondingly. 

We also performed a timing analysis of the ROSAT 
HRI observations. In both cases we used a 12 radius 
circular region to extract the source photons. This ra- 
dius encompassed 85% of the events in the ROSAT point 
spread function. The extracted data sets consist of 266 and 
521 counts, respectively. Photon arrival times were solar 
barycenter corrected, as implemented in the Extended Sci- 



entific Analysis Software System (Zimmerman et al. 1998). 

We analyzed the - 266 HRI (ROR 703348) photons 
from RBS 1223 for periodicities in the X-ray flux, but did 
not detect any significant signal using a Rayleigh Z\ test 
(Buccheri et al. 1983| ). However, a similar analysis of the 
second data set (ROR 704082, consisting of 521 photons) 
revealed a periodic modulation at 0.19394401 Hz (Fig. 
right panel). 

As Fig. shows, the significance of detection is not 
very high (Zf=23.2) owing to the small number of reg- 
istered photons. Moreover, the power spectrum is quite 
complicated; it consists of a central peak accompanied by 
many side lobes owing to the data gaps. This makes the 
computation of a credible region difficult, since the poste- 
rior PDF obtained by an application of the GL method is 
not unimodal, unlike in the case of the Chandra observa- 
tions. We determined the uncertainty of this period mea- 
surement of RBS 1223 through simulations, based on the 
measured pulsed fraction ~ 0.3 and count rate ~ 0.1s -1 . 
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Fig. 2. Periodogrammes of RBS 1223: Chandra ACIS-S (left) and ROSAT HRI (right) observations 
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Each of our 10000 simulated data sets (photon arrival 
times) was generated from a sinusoidal light curve and 
obeyed Poissonian statistics. Each simulated data set had 
the same characteristics as the observed data in terms of 
number of registered counts, spanned time, good time in- 
tervals, data gaps, and was analyzed exactly in the same 
way. The dispersion in the obtained period distribution 
was adopted as an estimate of the uncertainty in our mea- 
surement of the period of RBS 1223: SP = ±3.949 x 10" 5 s. 

An error estimate can be derived alternatively from a 
more general approach In general, the accuracy of a fre- 
quency determination in a periodogram depends on the 
total observation time T to tai, the number of measurements 



N and the signal-to-noise ratio (see e.g. Jaynes 1987, Brct- 
thorst |1988|, Kovacs |l98ll) . 



According to these authors, the theoretical accuracy 
for determination of the frequency of a single steady sinu- 
soid (di = A cos (2tt fti)+e(U), single sine wave plus noise 
e(ti) of standard deviation a) is 

5f= ^^Hz, 

' A ■ T tota i ■ y/N 

Hence, for the ROSAT observation of RBS 1223 (ROR 

704082) we get: T tota i ~ 70 ksec, N = 20 (number of bins 

in the phase folded light curve) and signal-to-noise ratio 

A/y/2a ~ 3 (see Fig. ||, right panel). This leads to an 

optimistic estimate of accuracy of the period of: 

dPoptimistic = ±4.4 X 10 5 S, 

while a pessimistic estimate is: 

8f = Hz (Rayleigh resolving power) 

2T tota / 

3 Ppessi?nistic — ±3.8 X 10 S. 

Thus, we finally derive the period of the detected pul- 
sations in the ROSAT HRI observations 
Pros at = 5.1561274 ± 3.8 x 10~ 4 s 

at the epoch of 50824.311057318 MJD. 

This determination allows to calculate a period deriva- 
tive and its uncertainty 
P = 1.35±g;|? x 10" 11 ss-\ 



if we adopt the pessimistic error estimate of the ac- 
curacy for period determination of RBS 1223 of ROSAT 
HRI observations (SP pessim istic)- 

2.2. Spectral analysis 

Since the Chandra X-ray observatory ACIS-S detectors 
have intrinsic energy resolution and the source spec- 
trum has sufficient number counts (count rate ~ 0.82 ± 
0.02 cts s _1 ), we were able to examine the spectral prop- 
erties of RBS 1223. After extracting (CIAO tool "dmex- 
tract") the spectrum of RBS 1223 and the nearby back- 
ground, we fitted the spectrum with an absorbed black- 
body model using the XSPEC VI 1.0 package. 

The ACIS-S instrumental response is not yet known 
in the energy range < 0.45 keV (see caveats page at 
http : II asc.harvard.edu/cda/ caveats. html). Hence, the 
absorption column density (Ah) cannot be well con- 
strained by the Chandra observation. The ROSAT PSPC 
spectra based on All-Sky Survey observation extends down 
to the 0.1 keV. We therefore in our fitting analysis first 
of all used only ROSAT PSPC data in order to determine 
Nn and then fixing that value to fit Chandra ACIS-S data 
in the energy range 0.45 1.5 keV, since no photons above 
1.5 keV were detected. 

The extracted spectrum and corresponding fit of 
ROSAT All-Sky Survey observation is shown in Fig. [l| 
The best fit parameters are a column density of Ah = 
(2.43 ± 1.06) x 10 20 cm" 2 , and a temperature kT BB = 
(113.7 ± 11.8) eV with formal 1 a errors. 

With the fixed value of Ah, the Chandra ACIS-S spec- 
trum of RBS 1223 can be represented by blackbody radia- 
tion with temperature kT BB = (90.6±0.8) eV with formal 
1 a errors (Fig. ||). 

The total galactic column density in the direction of 
RBS 1223 is A H , ga i = 2.1 x 10 20 cm" 2 , according to the 
Dickey & Lockman, ( 1990 ) maps of HI, averaged into 1° x 
1° bins (routine nh in FTOOLS). 
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Fig. 3. Phase folded light curves of RBS 1223: Chandra ACIS-S (left) and ROSAT HRI (right) observations 
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ROSAT All-Sky Survey spectrum of RBS 1223 
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Fig. 4. A blackbody fit and 68%, 90% and 99% confidence contours (Ax 2 = 3.5, 6.3, and 11.3, respectively) to the 
ROSAT PSPC spectrum of RBS 1223 



Chandra ACIS-S spectrum of RBS 1223 
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Fig. 5. A blackbody fit to the Chandra ACIS-S spectrum 
of RBS 1223. A best-fit blackbody temperature is kT = 
(91 ± l)eV (A H = (2.43 ± 1.06) x 10 20 cm" 2 ). 



In Fig. we show the "on-pulse" spectrum which is 
defined by averaging over the phase interval (0.2-0.6), i.e. 
it includes the maximum of the light curve in Fig. |3j. For 
comparison, we also show the "off-pulse" spectrum, ob- 



tained by averaging over the remaining phases. On- and 
off-pulse spectra were rebinned such that each spectral bin 
contains at least 25 photons. A spectral fit, again with 
fixed column density Ah = 2.43 x 10 20 cm~ 2 , yield best- 
fit blackbody temperatures of fcT(on) = 93 ± 3eV and 
fcT(off) = 88 ± 3eV with formal 3 a errors, indicating a 
trend that the source is hotter at maximum phase (Fig. ph. 



3. Discussion 

The detection of a periodic signal with P ~ 5.16 s proves 
that RBS 1223is a neutron star. Our determination of P ~ 
(0.7 - 2.0) x 10~ n ss" 1 shows that RBS 1223 is spinning- 
down very rapidly. Assuming a spin-down law of the type 
/ ~ — /" (/ = 1/P is the pulsar rotation frequency, / is 
its derivative, and n = f f /(f) 2 is the "braking index") 
with a constant magnetic moment gives the age of the 
pulsating neutron star, 



P 



(n- 1)P 

where Pq is a period of the pulsar at birth. 
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Fig. 6. The on/off pulse spectra of RBS 1223 obtained 
from Chandra ACIS-S observation. The on-pulse spectrum 
(upper) corresponds to the phases 0.2 to 0.6 on Fig. ||, 
while the off-pulse spectrum (lower) is from averaging over 
the remaining phases. A best-fit blackbody temperatures 
are kT (on) = 93 ± 3eV and kT(o//) = 88 ± 3eV (JVh = 
2.43 ± 1.06 x 10 20 cm- 2 ). 



Usually, an oblique rotating vacuum dipole m odel i s 
assumed, for which n = 3 (Manchester & Taylor 1977 ). 
If Pq is much smaller than P, the expression for the age 
determination reduces to r = Pj (2P) = r c , the character- 
istic age of a pulsar. 

Thus, the characteristic age of RBS 1223 is r c ~ (6000 
- 12000) yr. 

Assuming a neutron star radius of 10 6 cm and a mo- 
ment of inertia of 10 45 gcm 2 , the dipole magnetic field 
strength at the pole is 



B pole = 6.4 x 10 19 VPP ~ (3.5 - 6.5) x 10+ 14 G, 

well above the so-called "quantum critical field" value 

B c = = 4.4 x 10 13 G. This indicates that RBS 1223 
could be a highly magnetized young neutron star. 

There has been growing recognition in the literature 
of a population of highly magnetized neutron stars. The 
circumstantial evidence for this class comes from studies of 
soft gamma— ray repeaters (SGRs) and anomalous X-ray 
pulsars (AXPs). 

SGRs are neutron stars showing multiple bursts of 
gamma-rays with rather soft spectra (e.g. Hurley 200C for 
a recent review) . The spin periods of SGRs are clustered 
in the interval 5-8 s. They all appear to be associated with 
supernova remnants, which limits their average age to ap- 
proximately 20 kyr. The angular offsets of the SGRs from 
the apparent centers of their associated supernova rem- 
nant shells indicate that SGRs are endowed with space 
velocities > 500 km s -1 , greater than the space velocities 
of most radio pulsars (Cordes & Chernoff 199? ). 
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Fig. 7. Plot of P versus P for radio pulsars (dots), 
anomalous X-ray pulsars (AXPs), and soft gamma-ray re- 
peaters (SGRs). The dotted line shown between the lines 
for B = 10 13 and 10 14 G indicates a hypothesize d ap- 
proximate theoretical boundary (Baring & Harding 1998 ) 
separating radio-loud and radio-quiet neutron stars due 
to effects relating to magneti c field s close to the critical 
field B c (from Camilo et al. |2000[ ). RBS 1223 is identi- 
fied by filled box. A solid line is indicating where mag- 
netic field energy and rotational energy losses are equal 
(E mag = -±;B pole B pole R 3 and E rot = 4iT 2 Iff, Heyl & 



Kulkarni 199S) 



6 — 12 s, and have similar persistent X-ray luminosities 
as the SGRs (~ 10 35 ergs s" 1 , see e.g. Stella et al. |1998| ). 
Most of the AXPs appear to be associated with super- 
nova remnants, and therefore they are also thought to be 
young neutron stars like the SGRs. The spin periods of 
both AXPs and SGRs are increasing with time (spinning- 
down), and show no evidence for intervals of decreasing 
spin period (spin-up), although the spin-down rates of 
many of the SGRs and A XPs a ppear to be variable or 
"bumpy" (e.g. Woods et al. |2000| ). 

In addition, it is worth noting that some of the SGRs 
and AXPs are also located well outside apparently associ- 



AXPs are similar to SGRs in that they are radio quiet 
X-ray pulsars with spin periods clustering in the range 



ated s upern ova remnants (Marsden et al. 2001). Gaensler 
et al. ( 2001 ) argue that they may not be associated at all 
with these supernova remnants. 

Our values of P and P are suggestive, that RBS 1223 
might belong to either one of the two classes of unusual 
X-ray pulsars (Fig. 0). 

However, it is not excluded also that RBS 1223 may 
have a conventional (~ 10 12 G) magnetic field strength 
and identified as a neutron star in the propeller phase, 
i.e. powered by inflow of matter from a disk around it 
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(van Paradijs, Taam & van den He uvel |1995| , Chatterjee, 
Hernquist, & Narayan |2000| , Alpar |2001 ) 

Due to the lack of an identification at any other wave- 
length, the distance of RBS 1223 and hence its luminosity 
is quite uncertain. 

Assuming a size of the X-ray emitting region of 3-10 
km, our best fit spectral parameters indicate a distance of 
450-1500pc (~ 1.2 kpc assumed by Motch fl200C| )). 

Our determination of the absorption column density 
iV H = (0.6 - 6.6) x 10 20 cm- 2 (3er) (see Fig. |) poses no 
further constraint on the distance of RBS 1223. The ex- 
istence of diffuse high galactic latitude cirruses ( at the 
distance of 100-200 pc, see Magnani et al. |1996[ ) C9363 
and C9383, which extend to the position of RBS 1223 
(clearly seen in both IRAS 60/ik and IRAS 100/ik images, 
and having column density N-^ ~ (3 — 4) x 10 20 cm~ 3 , 
Reach et al. 1993| ) suggests that RBS 1223 might be a 
nearby (~ 100-200 pc) or rather distant (~ 700-1500 
pc) object. 

However, the absence of any obvious supernova rem- 



nant in the general direction of this source (Green 2000) 
and the fact that such a close and young supernova rem- 
nant would not remain undetected in the X— rays, suggests 
the possibility that the source has a very large transverse 
velocity and it may have migrated far from its original 
position. 

In this context it is worth mentioning that RBS 1223 
was not detected as a radio sourc e dur ing the VLA FIRST 
survey at 1.4 GHz (White et al. 1997 ) indicating a radio- 
quiet nature. The upper limit of radio flux at the source 
position is 0.94±0.138 mJy/beam. 

The position of RBS 1223 at the high galactic latitude 
of b 11 = +83° imposes an additional restriction on its dis- 
tance. If RBS 1223 would be at a distance ~ 1000 pc and 
would be born in the galactic plane, it must have traveled 
with a transverse velocity v t ~ 1000 km s _1 for about 
t ~ 10 6 yr. On the other hand, blackbody temperatures 
in excess of 50 eV imply ages younger than ~ 10 6 yr for 
standard cooling curves (see Tsuruta 199S ) . This suggests 
a much smaller distance than 1 kpc. 

Assuming a distance of the source of ~ 1000 pc and 
kTBB — 91eV, we determine the unabsorbed luminosity 
to be 

L x ~ 4.1 x 10 32 ergs s" 1 . 



4. Conclusions 

We analyzed Chandra ACIS-S and an archival ROSAT 
HRI observations of the isolated neutron star candidate 
RBS 1223 and found a period P ~ 5.16 s. We also found 
that the neutron star spins down at P ~ (0.7 — 2.0) x 
10~ n s s™ 1 with some remaining uncertainty on the spin 
down rate. Confirmation of our result requires further X- 
ray observations of RBS 1223. 

The nature of RBS 1223 remains uncertain, even if the 
spin down of the pulsar can be confirmed independently. 



It might be a young neutron star either with a very strong 
magnetic field (magnetar) or a neutron star with a more 
conventional magnetic field of ~ 10 12 G in the propeller 
phase powered by inflow of matter from a surrounding 
disk. 

The absence of a nearby supernova remnant and the 
'characteristic age' of a pulsar suggest a high space veloc- 
ity to bring RBS 1223 away from its progenitor supernova 
remnant. 

Additional constraints on its nature can be derived 
from detailed X-ray spectroscopy using upcoming XMM 
data and/or an X-ray proper motion study which we have 
already initiated. 
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